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Abstract. We review the current status of Andreev reflection spectroscopy on 
the heavy fermions, mostly focusing on the case of CeCoIns, a heavy-fermion 
superconductor with a critical temperature of 2.3 K. This is a well-established technique 
to investigate superconducting order parameters via measurements of the differential 
conductance from nanoscale metallic junctions. Andreev reflection is clearly observed 
in CeCoIus as in other heavy-fermion superconductors. Considering the large mismatch 
in Fermi velocities, this observation seemingly appears to disagree with the Blonder- 
Tinkham-Klapwijk (BTK) theory. The measured Andreev signal is highly reduced to 
the order of maximum ~ 13% compared to the theoretically predicted value (100%). 
The background conductance exhibits a systematic evolution in its asymmetry over 
a wide temperature range from above the heavy fermion coherence temperature 
down to well below the superconducting transition temperature. Analysis of the 
conductance spectra using the extended BTK model provides a qualitative measure 
for the superconducting order parameter symmetry, which is determined to be d^2_y2- 
wave in CeCoIns. It is found that existing models do not quantitatively account for 
the data, which we attribute to the intrinsic properties of the heavy fermions. A 
substantial body of experimental data and extensive theoretical analysis point to the 
existence of two fluid components in CeCoIns and other heavy-fermion compounds. A 
phenomenological model is proposed employing a Fano interference effect between two 
conductance channels in order to explain both the conductance asymmetry and the 
reduced Andreev signal. This model appears plausible not only because it provides 
good fits to the data but also because it is highly likely that the electrical conduction 
occurs via two channels, one into the heavy electron liquid and the other into the 
conduction electron continuum. Further experimental and theoretical investigations 
will shed new light on the mechanism of how the coherent heavy-electron liquid emerges 
out of the Kondo lattice, a prototypical strongly correlated electron system. Unresolved 
issues and future directions are also discussed. 
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1. Introduction 

It is well understood that the resistance minimum observed in metallic alloys diluted 
with small amount of magnetic impurities (e.g., Au-Fe or Cu-Fe) is caused by the spin- 
spin scattering between conduction electrons and local moments [1]. At temperatures 
approaching absolute zero, the magnetic moments of impurity atoms are screened by 
the conduction electron spins, prohibiting the divergence of the resistance due to the 
logarithmic scattering rate. This single impurity Kondo effect has been extensively 
studied in a variety of condensed matter systems [2|, [3] . If localized magnetic moments 
form a dense ordered array, i.e., Kondo lattice, the underlying physics becomes much 
more complicated and its understanding has been a great challenge to theorists and 
experimentalists. According to Doniach's Kondo lattice model [1], the fate of the local 
moments in the absence of direct dipolar interaction depends on the relative strength 
of two energy scales: one for Kondo interaction (/cb^k = -De"^/^'^^^^-'); and the other 
for Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction (/cb^rkky = J^N{Q))^ where 
k-Q is the Boltzmann constant, D is the bandwidth, J is the exchange coupling and 
A^(0) is the conduction electron density of states (DOS). The ground state will be an 
antiferromagnetically ordered one if Trkky > ^k, whereas a Fermi liquid of heavy 
electrons emerges in the opposite case. This heavy electron liquid, also called a Kondo 
liquid, does not exhibit spontaneous magnetism. Many theoretical and experimental 
efforts have been made towards a microscopic understanding of such phase diagrams 
[HI E] , particularly near quantum critical points [3, El E] . The microscopic mechanism 
for the formation of heavy electron bands remains unknown as the topic of the most 
fundamental importance to this field. A key question is how the localized f-electrons 
embedded in a quantum sea of conduction electrons acquire the itinerancy over the 
system. One scenario is that the f-electrons become delocalized via hybridization with 
the conduction electrons: a highly dynamical and collective process. This collective 
compensation of Kondo lattice spins by conduction electron spins forms new Bloch 
states and results into a narrow renormalized band of heavy electrons [6] . 

The majority of heavy electron materials are based on some rare earth (mostly 
Ce and Yb) or actinide (mostly U) elements that have valence electrons in the highly 
localized 4f or 5f orbitals. The overlap between neighboring f orbitals is negligible since 
they are highly confined to the nuclei over a length scale less than the interatomic 
spacing. Among dozens of Ce-based heavy-fermion compounds discovered so far, the 
1-1-5 family CeMIus (M = Co, Rh, Ir) have added to new avenues for the study of rich 
and novel physics including quantum criticality, the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) phase transition and the interplay between superconductivity and magnetism 
(see [Tot [TT] and references therein). The Ce 1-1-5 family is a subgroup of a larger 
class: CenMmIn3„+2m- Here, the Ce atoms occupy basal planes of a tetragonal crystal 
structure. Measurements of field-oscillatory magnetization (de Haas-van Alphen effect) 
in CeCoIns have identified multiple bands crossing the Fermi level [l2l [13]. Major 
Fermi surfaces consist of warped cylinders for heavy electron and hole bands and small 
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ellipsoidal pockets for light holes fTH IT3] . The resistivity maximum that signifies the 
beginning of heavy fermion coherence (T*) occurs around 45 K and a rapid drop is 
followed with decreasing temperature [H] . The emerging heavy electron liquid is found 
to disobey the Landau theory of Fermi liquid as evidenced by the T-linear (below ~ 20 
K) temperature dependence of the resistivity and the InT dependence of the electronic 
specific heat coefficient, to name a few [10]. This non- Fermi liquid behaviour implies 
proximity to a quantum critical point [HI |16l [T7] . Superconductivity in CeCoIus sets 
in below 2.3 K and the large specific heat jump at Tc indicates that heavy electrons 
participate in the Cooper pairing |14j . 

Most of known heavy-fermion superconductors exhibit unconventional pairing 
symmetry and gap structure [6]. For CeCoIus, the power law dependences of electronic 
and heat transport and the thermodynamic properties point to the existence of line 
nodes on the Fermi surface ^U\. Magnetic field-angle dependent thermal conductivity 
[T8] and specific heat [19] measurements exhibit four-fold oscillations in the basal plane, 
further supporting that the gap has d-wave symmetry [18], [IHl [20]. However, these 
results are seemingly contradictory for the precise node locations. Thermal conductivity 
indicates dx2_y2 [TBJ and specific heat supports d^y [I9]. There have been theoretical 
analyses which explain the controversy [21]. Not only is the interpretation of these 
experiments complex [21], but also they intrinsically cannot provide phase information 
of the order parameter. An unambiguous determination of the pairing state between 
dxy and d^2_y2 is of crucial importance to elucidating the pairing mechanism (e.g., see 

[22]). 

Tunneling spectroscopy as a probe of the superconducting order parameter has 
proven extremely useful and thus has been extensively adopted for both conventional 
[23] and unconventional [211 |25] superconductors. This technique was instrumental 
in confirming the predictions of the Bardeen-Cooper-Schrieffer (BCS) theory of 
superconductivity in conventional systems by measuring the quasiparticle DOS and 
phonon modes mediating the Cooper pairing ([23] and references therein). However, 
such a direct measurement of the order parameters in heavy-fermion superconductors 
has proven much more difficult compared to conventional superconductors and even to 
other unconventional superconductors, where planar tunneling and scanning tunneling 
spectroscopy (STS) have played key roles. This is due to their unique materials 
properties such as extremely low T^,, short coherence length, high reactivity, dependence 
of superconducting properties on the structural disorder, difficulty of thin film 
growth and non-cleavability. As an alternative technique, Andreev refiection (AR) 
measurements based on point-contact methods have been frequently adopted [26 | 127 ] l28] . 
This technique, called point-contact spectroscopy (PCS), has been successful to some 
extent but also has raised some issues including non-spectroscopic effects and how 
AR is possible at a normal- metal/heavy-fermion superconductor interface. The non- 
spectroscopic effects can originate from a contact being in the non-ballistic regime due 
to short mean free path, inter-grain Josephson coupling, etc, as will be discussed below. 

We review the current status of experimental and theoretical studies of AR in 
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heavy-fermion systems, focusing on CeCoIns. In section 2, basic theories on the charge 
transport in normal-metal/superconductor (N/S) junctions are reviewed. Basics of AR 
measurements using PCS technique are described in section 3 and technical issues in 
PCS are discussed in sections 4. Section 5 is devoted to the discussion whether an 
observation of a zero-bias conductance peak can be taken clS cl SI gnature of the sign 
change in the order parameter. In section 6, our experimental results are presented in 
terms of the superconducting order parameter symmetry. How AR is possible in heavy 
fermions is discussed in section 7. Possible origins of the conductance asymmetry, a 
common feature in all our measurements on the Ce 1-1-5 family, are addressed in section 
8. Our conductance modeling based on a possible Fano resonance in the Kondo lattice is 
described in section 9. Unresolved issues and future directions are presented in section 
10. Finally, conclusive remarks will follow. 

2. Charge transport across an N/S interface 

While trying to explain an unexpected increase of thermal resistance in the intermediate 
state of tj^e-I superconductors, in 1964 Andreev [29] discovered that there must be an 
additional scattering at interfaces between normal and superconducting domains. Its 
microscopic process was precisely derived by solving the boundary condition problem 
using Bogoliubov-de Gennes equations. The hidden scattering is retroreflection of an 
electron as a hole as it approaches an N/S interface from the N side. This particle- 
hole conversion process, called Andreev reflection, is essentially a scattering off a 
superconducting pair potential, during which energy, momentum and spin are conserved. 
The missing electronic charge is transferred as a Cooper pair into the superconductor, 
doubling the conductance within the energy gap, A. A schematic of this process 
is depicted in figure [H If the normal metal is replaced with a ferromagnet, AR is 
suppressed since the DOS for the opposite spin band is severely reduced (see figure [T](b) 
with unbalanced spin bands in N). This phenomenon has been extensively utilized to 
measure the spin polarization of ferromagnetic materials |30j . 

At about the same time as Andreev formulated the occurrence of AR, de Gennes 
and Saint- James |31j considered a similar boundary condition problem in an N/S bilayer 
system and discovered that the quasiparticle DOS is modified by the induced bound 
states. Although such AR-related effects had been known for a long time (e.g., [32]). 
AR was not considered as a direct measure of the superconducting energy gap until 
a theory came out in 1982. While studying the charge transport properties of N/S 
and S/N/S junctions. Blonder, Tinkham and Klapwijk (BTK) came up with an elegant 
theoretical formulation to explain various characteristics in a unified model [HHl [34] . 
They considered a boundary condition problem with a delta-function potential barrier at 
an N/S interface, H6{x), and the quasiparticle energy in the superconductor determined 
by the Bogoliubov-de Gennes equations. For an incoming electron impinging on the N/S 
interface from the N side, there are four possible trajectories: normal reflection, AR as 
a hole, transmission into the same branch and transmission with a branch crossing. In 
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Figure 1. A schematic of Andreev reflection process, (a) A real space picture. N and 
S denote a normal-metal and superconducting electrode, respectively, (b) A picture 
in the E (energy) vs. DOS space when the S is biased positively. The dashed line 
indicates the Fermi level. An electron incoming from the N side, impinging on the N/S 
interface, is retroreflected as a hole and get paired with another electron with opposite 
spin and momentum, transferring as a Cooper pair into the superconductor over a 
region ^ (superconducting coherence length) near the interface. Note that energy, 
momentum, spin and charge (including a Cooper pair) are conserved in this process. 



the BTK theory, the differential conductance across an N/S interface is given by the 
following formula. 

^{V) = 2N{0)ev^S dE^ii^^^ [1 + A{E) - B{E)] , (1) 

where A^(0) is the DOS at the Fermi level, e is the electronic charge, vp is the Fermi 
velocity, S is the junction area, / is the Fermi distribution function, A{E) is the AR 
probability and B{E) is the normal reflection probability. 

A remarkable aspect of this model is that it can reproduce characteristic 
conductance features from AR to tunneling as a function of single parameter, Z = 
H/hvp, the dimensionless barrier strength. This is demonstrated in figure [2] for an s- 
wave superconductor as calculated by ([T]). For a purely metallic junction {Z = 0), the 
subgap conductance is doubled with respect to the normal state value, which is due 
to 100% probability of AR for E < A. With increasing Z, a double-peak structure 
develops and the conductance at zero bias decreases, both features due to the energy 
dependence of AR probability. For a large Z value, say, Z > 3, AR is almost suppressed 
and the conductance shape is dominated by single particle tunneling. In this limit, the 
tunneling conductance maps out the quasiparticle DOS which, in the strong coupling 
case, may contain information on the bosonic modes mediating the Cooper pairing |35j . 
Since AR is a scattering off a pair potential, it is sensitive to the phase coherency of 
the Cooper pairs |36l|37j. Thus, AR measurements on unconventional superconductors 
may give different information than tunneling spectroscopy. This is an important issue 
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Figure 2. (a) Current vs. voltage and (b) normalized dl/dV vs. voltage characteristics 
calculated by the BTK formula for a normal-metal/s-wave superconductor junction at 
T=0. The dotted line in (a) is an asymptote to the Z—1.5 curve. = Rq{1 + Z"^), 
where Rq is the normal state resistance of a purely metallic junction. In the AR 
regime, an excess current is observed, which is the intercept obtained when the high- 
bias part of an I-V curve is extrapolated to the y-axis. Note that the BTK theory 
reproduces the transitional behaviour from AR to tunneling with a single parameter, 
Z, the dimensionless barrier strength. 

that has been debated actively in the case of high-Tj. cuprate superconductors [381 EH] • 
In an effort to analyze the conductance data obtained from Nb/Cu point-contact 
junctions, Blonder and Tinkham introduced a second term into the barrier strength 
parameter as follows [40] . 

Z,ff=[^' + (l-r)V4r]^/2^ (2) 

where r = Wp/ff, the ratio of Fermi velocities in the normal-metal (vp) and 
superconducting (ff) electrodes. They observed that the fitting parameter Zes ranged 
0.3 — 1 and interpreted that the minimum value of 0.3 is consistent with the Zes value 
calculated by ([2]) using r ^ 1.75 for Nb/Cu and assuming no insulating barrier {Z=0). 
Our experimental data on dozens of Nb/Au (r ^ 1.0 pX]) point-contact junctions over 
a range of sample-tip pressures have shown that Z^s > 0.2. This indicates that it 
may not be justified to extract the Fermi velocity mismatch contribution from the 
obtained Z^g values, contrary to their interpretation and to the claim by Deutscher 
and Nozieres [42] that one can determine the mass enhancement factor by comparing 
velocities obtained from PCS with known values. Regarding this issue, Lukic considered 
generalized boundary conditions in the BTK problem and argued that the effective mass 
is just a parameterization of the unknown microscopic parameters [13] . 

In heavy-fermion materials, the electronic mass is highly enhanced (by an order 
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of ~ 10 — 10'^) with a correspondingly reduced Fermi velocity. According to ([2]), a 
junction with a heavy fermion as one electrode and a conventional metal as the other is 
inherently in the tunneling regime (e.g., > 5 for Au/CeCoIns) and thus AR cannot 
occur. However, AR has been frequently observed in such junctions as we present our 
data below and as observed by others ([27] and references therein). Deutscher and 
Nozieres addressed this discrepancy by proposing that the boundary conditions are not 
affected by the mass enhancement factor [12]. These issues are further discussed in 
section 7. 

3. Basics of point-contact spectroscopy 

Detailed reviews on this topic have been published, e.g., [28]. Here we only present 
fundamental concepts and elements of PCS. Let us begin by considering electrical 
conduction across a metallic contact between two identical metal electrodes. Assume 
the two metals are connected only through a constriction of radius a and otherwise 
electrically isolated. If a bias voltage is applied across the junction, the resulting 
conductance will depend on how the electric field distributes in the contact region. 
Here, the most relevant length scales are the electronic mean free paths, elastic (/) and 
inelastic (/in), in comparison to the contact radius a. If 2a ^ electrons gain or 
lose energy only when crossing the interface and their velocity change is proportional to 
the applied voltage. This gives rise to an ohmic current- volt age characteristic and the 
contact resistance is given by the following formula, known as Sharvin resistance |44j : 



where p is the resistivity of the metal. This limit (2a ^ Z, /in) is called ballistic or Sharvin 
limit. Note that the resistance depends only on the contact area since the product pi is 
constant in the simplest Drude picture. In the opposite limit, namely, if 2a ^ I, /in, the 
electron distribution would be that as in the bulk. Then, the point-contact resistance 
is given by the formula for a bulk sample of length 2a and diameter 2d [35] , except for 
a geometric factor: 



This limit (2a ^ /, /in) is called the thermal or Maxwell limit and a local Joule heating 
can predominate the electrical transport process. 

By solving Boltzmann transport equations for a point-contact junction with an 
arbitrary value for K = l/a, Wexler derived a generalized formula for the point-contact 
resistance as follows [46j. 




(4) 
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undergo multiple elastic scattering across the constriction. In the case of an N/S point- 
contact junction, this may reduce the AR signal. In figure [3], the three regimes are 
indicated in relation to the contact size and mean free paths. An extreme limit, not 
relevant to the PCS discussed in this paper, can occur if the contact size becomes 
comparable to the Fermi wavelength, in which case quantum size effects dominate 
[171 EH]- This regime is reached in atomic size quantum point contacts |48j . 

If a point-contact junction is in the ballistic limit, energy-dependent quasiparticle 
scatterings appear as nonlinearities in the current-voltage characteristics, analysis of 
which gives spectroscopic information of the material under study. What makes this 
PCS possible is that the contact area remains cold due to a large mean free path 
compared with the contact size and any energy dissipation occurs away from the 
junction. In contrast, in the thermal regime, energy dissipation occurs within the 
contact volume, causing local heating. This PCS technique was first demonstrated by 
Yanson in 1974 in his second harmonic measurements of micro-shorted thin film tunnel 
junctions [19], from which he obtained phonon spectrum in Pb, similar to that obtained 
in the tunneling measurements by McMillan and Rowell [35]. Shortly after, Jansen 
and coworkers developed a more controllable method to make point-contact junctions 
[50] . This technique, called the 'spear-anvil' method, has been adopted most widely 
for PCS because of its simplicity. During its early days, PCS was mostly utilized to 
investigate phonon spectra in metals and alloys [51],|52]. Since the BTK theory came out 
[33] to explain the experimental data nicely [Ml HQ], PCS has been frequently adopted 
as an alternative tool to planar tunneling or STS for the study of superconducting gap 
structures of many novel and unconventional superconductors owing to its versatility and 
technical simplicity [27l [28], [53] . In the majority of PCS experiments on superconductors, 
the differential conductance vs. bias voltage data have been analyzed successfully using 
the BTK model to determine their superconducting order parameters. 

There exist many reports in the literature where the PCS data may not necessarily 
be spectroscopic. As discussed above, for purely spectroscopic measurements, the 
junction must be in the ballistic regime. A first check of this is to estimate junction 
size using ([5]) and compare with known values for the mean free paths . Because 
in general point-contact junctions are made by bringing two metals into mechanical 
contact on the nanoscale, data can also be affected by factors other than the nominal 
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Figure 3. Schematic of the PCS regimes, a is the radius of a point contact and I 
(Zin) is the elastic (inelastic) mean free path of electrons in the electrodes. Ballistic: 
2a < Diffusive: ^ < 2a < Thermal: 2a > l,kn- 
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contact size, such as contact pressure and geometry [SS], meaning that the nominal 
contact size being in the balhstic hmit alone cannot guarantee the spectroscopic nature 
of the data. The most important diagnostic is the reproducibility and consistency 
of the conductance spectra point to point, sample to sample and along different 
crystallographic orientations, as we discuss in the next three sections. 

4. Experimental methods 

Like for other class of superconductors, the spear-anvil method has been employed 
widely for the PCS investigation of heavy-fermion superconductors. Their much lower 
Tc's (typically less than 1-2 K) make its implementation complicated because of the 
inevitable involvement of a ^He or dilution refrigerator. In these purely 
mechanical approach of the tip to the sample based on the differential screw mechanism, 
which is widely adopted for ^He systems, is not suitable. Instead, a combination of 
mechanical and piezoelectric adjustments is more common, where a coarse approach is 
made by a fine-pitched screw and a fine adjustment is made by driving piezoelectric 
elements [56]. Although mechanically cut wires (typically Au or Pt-Ir) have been used, 
electro chemically polished tips provide better control over the tip size and surface 
morphology. This is an important technical issue since it could affect the contact 
geometry and thus the conductance spectra. 

It is also crucial to prepare smooth and clean sample surfaces in order to obtain 
spectroscopic data. Almost all materials investigated by PCS have surface layers that 
are not the same as in bulk. These layers could be precipitates from growth, surface 
reconstruction or due to oxide growth when exposed to air. In some cases, these layers 
can be penetrated by the metal tip during the tip approach, enabling one to probe bulk 
characteristics. Otherwise, the surface layer can act as a potential barrier, proximity or 
degraded layer. If the surface layer is thin, the potential barrier case will still provide 
spectroscopic data which can be analyzed by the BTK theory, but approaching, or in, 
the tunneling limit. In the other two cases, obtaining spectroscopic data is more difficult, 
particularly if the surface layer is thick on the order of the coherence length or longer. 

It is ideal if a superconductor under study can be prepared such that it has fiat 
and smooth surfaces along all major crystallographic orientations. Most unconventional 
superconductors are first produced in bulk (poly- or single-crystalline) forms. PCS 
on polycrystals can reveal some information on the gap values but the momentum 
direction is not identified or defined well, meaning that data analysis and interpretation 
might not provide crystallographic orientation-dependent information. In addition, 
because of a granular structure, measured data are vulnerable to artificial effects such as 
intergrain Josephson coupling [57]. Thus, single crystalline samples are highly desirable. 
A drawback, particularly in anisotropic materials, is that crystals are often grown as 
platelets only along the preferred surface orientations. Thin film samples can be grown 
along desired directions by adjusting growth conditions but heavy-fermion compounds 
are known to be extremely difficult to synthesize in thin film form. For PCS, we rely on 
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the high-quality single crystals that have been grown and the preferred as-grown surface 
of a single crystal is first investigated. For measurements along the other directions, 
broken edges are frequently used. Preparation of surfaces by polishing is rarely adopted 
due to the possibility of surface degradation during the polishing process. However, this 
approach for sample preparation has been successful for the AR study of CeCoIns |58] . 

The heavy-fermion compound CeCoIus has a tetragonal crystal structure with the 
Ce atoms occupying the basal plane. Many experimental reports have pointed to a 
d-wave symmetry of the Cooper pairing state [IDl HI UHl EO] • Considering its known 
Fermi surfaces [121 [13], the gap nodes are supposed to form lines along the c-axis, 
similarly to the case of cuprate superconductors. From the group theoretical point of 
view on the possible pairing state, there is a correspondence between real and momentum 
space regarding the gap node directions in d-wave superconductors [59]. Therefore, it 
is essential to make PCS measurements into the ab plane to probe along nodal and 
antinodal directions. More specifically, measurements into (100) and (110) surfaces 
should be made to identify the location of nodes. The largest face (typically ranges 
a few mm) of an as-grown CeColus crystal is perpendicular to the c-axis. Thus, this 
surface is naturally chosen for PCS along the c-axis. Since the thickness of single crystals 
ranges only a few hundred /im's, it is necessary to devise a method to hold a sample 
securely such that its in-plane direction is aligned with the tip axis. This can be achieved 
by embedding a crystal into a mold made of low-temperature epoxy and cutting and 
polishing it [58j. 

CeCoIn5 crystals are grown by the flux method using excess indium flux [H]. 
Because of the nature of this method, samples can contain indium precipitates, which 
may affect the PCS measurements. To eliminate this possibility, one slightly etches 
crystals using hydrochloric acid, noting that CeCoIus is also dissolved although at much 
slower rate than indium. Microscopically rough surfaces can therefore be produced, 
causing spurious effects. Again, to ensure the spectroscopic nature of the conductance 
spectra, reproducibility is the most important diagnostic. In particular, data should 
be compared along all major crystallographic orientations. In general, considering the 
extreme cleanness of CeCoIus (the electronic mean free path ranges several ~ /im's at 
low temperature [60]), it is highly feasible to perform PCS measurements in the ballistic 
regime. 

5. Zero-bias conductance peak vs. Andreev bound states 

Three independent groups have reported results of PCS on CeCoIus. Coll et al. [61] 
first reported two different types of conductance spectra obtained from point contacts 
along the c-axis: one with double-peak and zero-bias dip structure, and the other with 
a zero-bias peak and surrounding dips. The latter feature was interpreted as hinting at 
an unconventional symmetry (rf-wave was implied) . Rourke et al. [62] also reported two 
different types of spectra from nominally c-axis contacts: one with a zero-bias peak and 
hump structure, and the other with a double hump structure. These authors claimed 
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both (i-wave symmetry and multiple order parameter components. The validity of the 
claim for (i-wave symmetry based simply on the observation of a zero-bias conductance 
peak (ZBCP) has been questioned [55l [63] because many other origins for a ZBCP are 
also known. 

Figure H] shows normalized dl/dV curves calculated for a normal- met al/ci- wave 
superconductor junction using the kernel developed by Tanaka and Kashiwaya et al. 
[HH [65], |66] as a function of Zeg, the effective barrier strength in ([2]). Here, it is assumed 
that the injected electrons have a momentum distribution over the whole range of angle 
(90 deg.) around the junction normal. In practice, in the tunneling limit, tunneling 
electrons may have a much narrower momentum distribution, known as tunneling cone 
effect [67]. In the limit of a 90 degree tunneling cone, there is no difference between the 
c-axis, nodal or antinodal directions as there is an integration over the full half space 
of momentum, as shown in figures [H For an antinodal or c-axis junction, the cusp-like 
feature and linear slope near zero bias for purely metallic (Zgfj = 0) and tunnel junctions 
(Zeff = 5), respectively, exhibit the characteristic DOS of a ci-wave superconductor. This 
is just a qualitative difference from the s-wave case, where a flat conductance shape is 
seen due to a fully isotropic opening of the gap. As is increased from zero, a double- 
peak and zero-bias dip structure develops, which is also similar to the s-wave case (see 
figure [2]^b)). An interesting case is the nodal direction with a finite Zgff. The conductance 
shape varies in a strikingly different way, as displayed in figure [H^b). Beginning with the 
same shape for Zeg = as for the antinodal junction, the conductance curve becomes 
narrower and sharper with increasing Zeg, forming a ZBCP instead of a dip. This 
ZBCP originates from Andreev bound states (ABS) formed on the nodal surface due 
to a sign change of the ci-wave order parameter around the Fermi surface [68]. The 
ZBCP in tunnel junctions on hole-doped cuprate superconductors, well known to be d- 
wave, has been frequently observed [53], [69] and well understood theoretically [70]. The 
ABS occurs due to the constructive interference between an incoming electron and an 
Andreev-refiected hole and can be formed at zero energy on a nodal surface of a ci-wave 
superconductor because of the additional phase shift of vr [HE] [71] . The same shape for 
a nodal junction as for an antinodal junction when Z^q = can be understood as due 
to smearing of ABS in the absence of barrier [65l [66] . 

It has been frequently claimed that a ZBCP observed from PCS is a signature 
of ABS, thus, evidence for a sign change of the superconducting order parameter. 
However, as pointed out by us [55j, and Sheet and Raychaudhuri [63], just an 
observation of a ZBCP does not corroborate such claims. A ZBCP can also arise from 
other physical mechanisms without requiring a sign change such as local heating [72], 
intergrain Josephson coupling [57J, etc. Therefore, it is essential to perform diagnostic 
measurements and analyses in order to prove that a ZBCP is a real signature of ABS. In 
the high-Tc cuprate tunnel junctions, it is well established that ABS-originated ZBCPs 
are split either spontaneously or by applied magnetic field [53] [69] . One of the broadly 
accepted explanations is the Doppler shift of the bound state energy due to coupling to 
the superfiuid momentum [70]. The amplitude of an ABS-originated ZBCP is predicted 
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Figure 4. Normalized dl/dV vs. voltage characteristics calculated by the extended 
BTK formula for a normal-metal/d-wave superconductor junction at T=0. (a) The 
junction normal is along the antinodal direction. The conductance behaviour is 
qualitatively similar to the s-wave case except for the cusp-like (AR) and V-shape 
(tunneling) features near zero bias, which are due to a d-wave superconducting DOS. 
(b) The junction normal is along the nodal direction. Completely different behaviours 
are observed due to the surface Andreev bound states. 

to follow 1/T dependence, in contrast with the InT behaviour of Kondo scattering 
effect in tunnel junctions. In the case of metallic junctions as in PCS, it is proposed 
that the critical magnetic field over which the ZBCP can be split is much stronger 
|73j . Another important diagnostic measurement, particularly necessary for metallic 
junctions because the splitting may not be observed for the reasons above is to compare 
conductance spectra along different crystallographic directions. In the case of CeCoIus, 
given that c/-wave symmetry is most likely, measurements along both (100) and (110) 
directions are essential to differentiate between dx^-y^- and d^y-^wave symmetry. 

Neither of the reports by Goll et al. [61] and by Rourke et al. provided such 
diagnostic measurements, so origins of their observed ZBCPs remain unclear. Rourke et 
a/.'s claim on the existence of multiple order parameters relies on their observation 
of multiple hump structures in the differential conductance data. However, such 
multiple hump or peak and dip structures have been sometimes observed in PCS study 
of a known single gap superconductor. To prove that the hump positions can be 
taken as gap values for multiple order parameters, diagnostic measurements including 
crystallographic orientation and temperature dependences should be made along with 
reproducibility. In general, considering the band structure and Fermi surface topology 
[121 [13], it is possible that the superconductivity in CeCoIus has multiband nature, but 
it is a different issue to prove that measured PCS data actually reflect such features. 
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6. Spectroscopic evidence for dx2_y2 symmetry 

We have reported differential conductance spectra from Au/CeCoIns point-contact 
junctions along all three crystallographic directions [581 EH ES]. Figure [5](a) displays 
temperature evolution of the conductance along the c-axis over a wide temperature 
range. At high temperatures, the conductance curves are symmetric and flat, which is 
characteristic of simple metallic junctions over a small bias region. With decreasing 
temperature, an asymmetry in the dl/dV curves is developed. This conductance 
asymmetry appears to begin below the heavy fermion coherence temperature T* (~ 
45 K) f76] and increases with decreasing temperature down to Tc (2.3 K), below 
which it remains constant, as shown in figure [5](a). As is crossed, the conductance 
near zero-bias begins to be enhanced and this enhancement increases with decreasing 
temperature. In the literature, it has been frequently adopted to analyze conductance 
data by normalizing them with respect to those just above T^. Noting that our data 
below Tc exhibit nearly the same background shape, we normalized them against the 
data taken at 2.6 K. The result is plotted in figure [5](b). The temperature and bias 
voltage dependences and the shape of the conductance curves (BTK-like) near zero bias 
indicate the enhancement is due to AR. However, the magnitude of the AR conductance 
(~ 13%) is greatly suppressed in comparison to the theoretical prediction (100%) [33] . 

Our attempt to fit the normalized conductance spectra on the (001) CeCoIns 
junction using the extended BTK model by Tanaka and Kashiwaya et al. [Ml ESI [66] 
revealed that this model could not fully account for the data [71] . The best-fit curves 
deviate from the data substantially (figure [St^a)) and the temperature dependence of 
the fitting parameter F, smearing factor [77J, basically tracks that of the gap, which 
is unphysical (figure [6](b)) [71]. The large deviation around gap edges implies that 
the failure is not solely due to the large suppression of AR but also due to a large 
shift of the spectral weight [71]. This situation remains unchanged no matter which 
modified BTK model we adopt for s- or d-wave symmetry considering the Fermi surface 
mismatch in the two electrodes and the possible breakdown of Andreev approximation 
[ZU [THl [TU [80l [m [82]. All these results point to a necessity to develop a new 
theoretical framework to understand how AR is possible at all and why it is observed 
so reduced in heavy- fermion/conventional- metal interfaces. Even without the full 
theoretical framework, PCS can be used as a spectroscopic probe of the order parameter 
symmetry as we discuss below. 

We now address the superconducting order parameter symmetry in CeCoIus. As 
seen in the simulation in figure HI a cusp-like shape is predicted for d-wave symmetry if 
there is no barrier (Zeg = 0). The zero effective barrier strength condition is difficult, if 
not impossible, to achieve in conventional-metal/heavy- fermion-superconductor point- 
contact junctions because of some existing surface layer or even some effect of the 
Fermi velocity mismatch. We therefore assume Z^q is always finite. In this case, a 
cusp-like or peaked structure, even for small finite Zefr, would never be observed in the 
antinodal direction. There always exist smearing effects due to finite temperature, finite 
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Figure 5. (a) Differential conductance spectra of a (001) CeCoIns point-contact 
junction. Note the systematic evolution of the background conductance asymmetry 
and the enhancement of conductance near zero bias due to Andreev reflection, (b) 
Conductance curves normalized by the 2.6 K data. There is no dip in the raw data in 
(a), so the shallow dips around ± 1 mV at 400 mK in (b) are due to the normalization, 
not caused by local heating. The magnitude of AR conductance is on the order of ~ 
13%. After [TJ. 

quasiparticle lifetime and depairing. Thus, the interpretation of a flat conductance shape 
as we observed in (001) and (100) CeCoIns junctions is ambiguous even between s- and 
d-wave symmetry. This was reported earlier for (001) junctions [71], where a li-wave 
fit gives slightly better results. Therefore, it is essential to obtain conductance data in 
the afe-plane to determine the pairing symmetry. Figures [7](a) & (b) compare such data 
along the (100) and (110) directions. While both spectra exhibit similar background 
asymmetry as seen in the (001) junction, their detailed shapes in the sub-gap region 
are quite different from each other. Namely, the (100) data look flat, similarly to the 
(001) data, whereas the (110) data are cusp-like. A comparison of these experimental 
observations with theoretically predicted behaviours gives strong evidence for d-wave 
symmetry, in fact d^2_y2 [75j. Figures [Tt^c) & (d) show calculated curves in the small- Zeg 
limit using the d-wave BTK model [65]. The calculated conductance curves at Zcs = 
are identical, as discussed earlier. For finite Z^g, they develop completely different 
features: double-peak and dip for the antinodal junction vs. narrower and shaper peak 
at zero bias for the nodal junction. We conclude that the (100) data are consistent 
with antinodal junction behaviour with Zcs ~ 0.3 and the (110) data represent nodal 
junction behaviour with a smaller Zcs- Therefore, the superconducting order parameter 
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Figure 6. (a) Normalized dl/dV (□) of a (001) CeCoIns junction in figure[5Ub) and 

best fit curves ( ) by the s-wave BTK model, (b) Energy gap (A) and quasiparticle 

lifetime broadening factor (F) obtained from the s-wave fit. A predicted by the BCS 

theory ( ) and F ( ) obtained from the zero-bias conductance fit. Note that 

the temperature dependence of F is unphysical. After |74j . 

symmetry in CeCoIns is (i2.2_y2-wave, not (i^^y-wave. 

This spectroscopic result resolves the controversy of the node locations [T8l [T9l [2T] . 
This definitive determination of the pairing symmetry narrows down possible candidates 
for the bosonic mode involved in the microscopic pairing mechanism [22j. It is notable 
that AR measurements along different crystallographic directions provide spectroscopic 
evidence for the order parameter symmetry. This is possible via detecting the sign 
change albeit ABS are smeared to finite energy due to the higher junction transparency 
compared to tunnel junctions. Therefore, like ABS tunneling spectroscopy, PCS can 
provide phase-sensitive information in unconventional superconductors [831 El] • As for 
the anisotropy in the order parameter amplitude, the conductance widths are similar 
for the (100) and (110) junctions probably due to the large tunneling cone effect. This 
is an issue to be addressed more rigorously in the future. A cusp-like feature similar to 
that in our data is observed in GoU et a/.'s measurements into the afe-plane of CeCoIns 
but the precise crystallographic orientation was not reported [85] . 

As discussed previously, ABS are known to split spontaneously or by applied 
magnetic field. This splitting occurs due to the Doppler shift effect in the bound 
state energy arising from a coupling of the quasiparticle velocity {v-p) to the superfluid 
momentum (Ps), "^f ■ -Ps- According to Tanaka et al. [73], the critical field for splitting 
of ABS is proportional to the junction transparency. Point-contact junctions generally 
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Figure 7. (a) and (b) are normalized dl/dV of the (100) and (110) CeCoIns junctions, 
respectively, (c) and (d) are normalized dl/dV curves calculated by the d-wave BTK 
formula (r=0 and T=0) for antinodal and nodal junctions, respectively, (e) and (f) 
are magnetic field (parallel to the c-axis) dependence for the (100) junction at 400 mK 
and the (110) junction at 420 mK, respectively. After [75] , 



have high transparency, so sphtting of ABS would only be accomplished at very high 
field. This is completely consistent with our field-dependent measurements, as shown 
in figures [Tl^e) and (f). Other possibilities for the non-splitting of ABS have also been 
suggested, including the tunneling cone effect [HE] and atomic scale disorder [87] in 
cuprate tunnel junctions. 
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7. How Andreev reflection can occur in heavy fermions 

Our conductance data on CeCoIus is a clear example for AR in heavy fermions. 
Also, there have been many reports on the observation of AR in other heavy-fermion 
superconductors [26l[27l|88l|89l[90l[9ll|92l|93]. As mentioned previously, this 
experimental observation is in conflict with what the BTK theory predicts as in ([2]): 
AR should be completely suppressed due to the large mismatch in the Fermi velocities 
of quasiparticles. Thus, how AR is possible in heavy fermions is a longstanding issue. 
Deutscher and Nozieres |12] addressed this issue and proposed a scenario that the 
velocity (vp) of injected or extracted particles are not the same as the quasiparticle 
velocity (vp) as seen in specific heat or coherence length measurements. Namely, 

vp = vp ■ z, (6) 

where z is the mass renormalization factor and 1/2; = 1 + Ais the mass enhancement 
factor [94]. In heavy fermions, z is very small. Lukic considered generalized boundary 
conditions in the BTK problem and proposed that one should not interpret the 
parameter Z^s as representing materials properties but just as parameterization of 
the unknown microscopic details of the interface |43j. Even though these arguments 
provide a way to explain the discrepancy between the BTK theory and experimental 
observations, they do not account for why AR is strongly reduced in heavy fermions. It 
is our conjecture that this issue might be intimately connected to why the conductance 
shows an asymmetry, so we focus on this issue in the next two sections. 

8. Conductance asymmetry 

In our PCS data on CeCoIns, there are two prominent and common features that are 
not seen in other class of superconductors: highly suppressed AR and background 
conductance asymmetry. It is noted that these features have also been observed in 
other reports on the PCS of heavy-fermion superconductors [271 ESI ES, 196] . 

Although our qualitative analysis of the point-contact conductance spectra using 
the d-wave BTK model gives spectroscopic evidence for the superconducting order 
parameter symmetry in CeCoIus, it does not give us precise information on the gap 
amplitude. As discussed previously, existing models taking the Fermi surface mismatch 
and possible breakdown of Andreev approximation in heavy fermions into consideration 
do not quantitatively account for our data. This is an indication that a crucial element 
is missing in the currently available models. This situation is in striking contrast to 
those cases in the other class of superconductors, where the BTK models have been 
so successfully applied |2Hl [53]. We attribute the origin of this failure to the intrinsic 
properties of heavy fermions, more specifically, the coexistence of multiple electronic 
components [76l [97] [98]. 

In figure M^a), it is found that the temperature dependence of the conductance 
asymmetry qualitatively follows that of the spectral weight for the coherent heavy 
electron liquid reported by Nakatsuji, Pines and Fisk [76]. Also, there is a systematic 
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Figure 8. (a) Temperature dependence of the background conductance asymmetry, 
defined as tiie ratio between conductance values at -2 mV and at +2 mV from the 
(OOf) CeCoIus data in figure [5] T* is the heavy fermion coherence temperature. The 
inset is a semi-logarithmic plot of the same data, (b) Temperature-dependent evolution 
of the normal state conductance. After [75l. 



evolution of this background conductance peak position as shown in figurel^b): it moves 
from the far-negative bias side towards zero bias with decreasing temperature. Yang and 
Pines [98] recently reported a universal scaling behavior between various experimental 
data including our conductance asymmetry in figure [8](a) and the spectral weight for 
the emergent heavy-electron liquid. 

According to Harrison's theorem [99], the electronic DOS of a simple metal cannot 
be measured using tunnel junctions since it cancels out the velocity factor in the 
conductance kernel. This theorem appears to hold in many cases where the single 
particle picture is valid. Applying this argument to point-contact junctions consisting 
of simple metals, the electronic DOS is not expected to be measured from such 
measurements. The conductance of a metallic contact is given by 

^(-'-//-'-)^^--' 

where v is the velocity, N{E) is the electronic DOS of the counter-electrode, / is the 
Fermi distribution function and dQ is the differential solid angle [27]. Over a small 
bias region, the conductance is expected to be flat. Beyond that region, it is expected 
and observed to be curved downward due to increased scattering, say, with phonons, as 
the energy increases [47], and this shape is not attributed to the DOS. The downward 
curvature in metallic junctions is contrary to the upward curvature in tunnel junctions 
|47] . Now, the question is whether Harrison's theorem still holds in the case where the 
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many-body interaction is not negligible, as in heavy fermions. 

The observed conductance asymmetry in our data on CeCoIus might be strongly 
tied to this issue. As discussed in detail in |100j . the conductance asymmetry cannot be 
explained by the models proposed in the literature. Certainly the asymmetry indicates 
it is easier to remove electrons from than to add to the heavy-fermion electrode CeCoIns. 
However, a simple argument that it is more difficult to add than to remove an electron 
from the f-orbital does not work because the f-electrons do not remain localized, but are 
itinerant, in heavy-fermion compounds. The argument that competing orders cause a 
conductance asymmetry |101j doesn't appear relevant to our case because it is observed 
in all three members of the Ce-based 1-1-5 heavy fermions |100] . no matter how far they 
are from the competing region. The persistency of conductance asymmetry in magnetic 
field up to 9 T (see figures [7](e) & (f)) rules out a connection to the argument based on 
non- Fermi liquid behaviours [I02j . In the past, the conductance asymmetry observed 
in PCS of heavy fermions was frequently attributed to the large Seebeck coefficients of 
heavy fermions [1031 11041 1105j combined with contact's being in the thermal regime. 
However, this model doesn't agree with our experimental data on 1-1-5 heavy fermions 
|100] . A simple counter argument to this scenario is why the asymmetry persists in the 
superconducting state as seen in figure [5](a), where the Seebeck effect disappears [106] 
and, thus, the asymmetry should, too. 

Anders and Gloos proposed that both the conductance asymmetry and reduced AR 
in heavy-fermion superconductors could be caused by the strongly energy-dependent 
quasiparticle scattering and DOS of heavy fermions [107] . This theory does account 
for our observed conductance features but, since the calculations are done in a Green's 
function formalism, the microscopic physical mechanism is not obvious. Nowack and 
Klug [108] used Boltzmann transport equations to deal with electronic scattering by an 
energy-dependent DOS of heavy fermions. In this model, the asymmetry arises from 
the DOS centered at a finite energy. We note that the DOS effect of heavy fermions is 
invoked in both models. 

9. Fano interference effect and conductance modeling 

We have proposed |75] a two-channel conductance model based on the two-fluid picture 
[76] and originally assuming a Lorentzian form of the heavy electron DOS. Here, the 
two channels are into the superconducting heavy electron liquid and a possible normal 
conducting light electron liquid [109] . It was assumed that the two channels are 
independent of each other. In theory it was argued that superconductivity should be 
induced to the light normal electron liquid due to the proximity effect [110] but this issue 
is not yet settled. Our experimental data were found to accurately flt to this model at 
lowest temperatures but the data couldn't be flt with the model at higher temperatures; 
the flt failing more severely with increasing temperature. This is technically because the 
background conductance shape is not strictly Lorentzian. Instead, it is observed [25j to 
resemble a Fano line shape [111] , implying that the two conductance channels need to 
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be entangled instead of being independent. Motivated by our two-channel conductance 
model, Araiijo and Sacramento recently formulated a BTK model adopting the two- fluid 
picture for heavy fermions [112|. In this model, AR is reduced essentially due to the 
existence of a normal conduction channel, similarly to our model, but they claim that 
the two channels should be allowed to mix via an interference term instead of being 
independent. 

The Fano resonance originally discovered in electron-helium inelastic 

scattering cross sections, has been frequently observed in a variety of condensed 
matter and other physical systems. In essence, it is a manifestation of an interference 
effect between waves that have passed through two different paths containing discrete 
and continuum states, respectively. This interference causes an asymmetry in the 
resonance spectra whose line shapes can be reproduced by the Fano formula [lllj . 
Recently, this phenomenon has been a subject of intensive investigations in quantum 
dot experiments [3l [ml [IH [Iia [M [HTl [IH] and STS of single adatoms and 
molecules [Il9l [EOl [Ml [1221 [1231 [IMl [ESI [ESI [EH [EHl [M [I30l [ml [I32l [188] . 
For example, Kobayashi et al. observed systematic evolution of the asymmetry in the 
differential conductance across a mesoscopic ring containing a quantum dot in one 
arm and a continuum in the other arm |1171 I118j . Madhavan et al. investigated 
conductance behaviours of single adatoms deposited onto metallic substrates using 
STM |119l 1124] . The measured differential conductance showed featureless shapes from 
non-magnetic atoms, whereas highly asymmetric curves from magnetic atoms. Their 
proposed conductance formula including the Fano resonance term successfully accounts 
for the measured data |124j . Here, the interference occurs between tunneling electrons 
into the conduction band and into the local d orbital of an adatom. 

The original Fano formula is given as: 

m - (8) 

where q-p is the Fano parameter, e = {E — Eq) / 1- and Eq and A are the resonance energy 
and width, respectively. Using this formula, Fano lines are simulated as a function 
of qp, as plotted in figure [9](a). In the extreme limit of Igpl —>■ oo, the line becomes 
symmetric and represents a Lorentzian resonance. In the opposite limit, i.e., \qp\ —>■ 0, 
it corresponds to an antiresonance. In between these extremes, the line is asymmetric 
with a peak positioned at e = l/qp. 

Adopting a similar formula to that used by Madhavan et al. |124j , we were able to 
fit the data and extract useful information [100] . We set up a model that the differential 
conductance is determined by the following equation containing the Fano formula ([8]) 
as the kernel. 

where C is a prefactor for the Fano resonance contribution and Gq is a constant 
conductance. An example of the fit using this formula is shown in figure [H Here, 
the used fit parameters are: gp = -2.14, Eq = 2.28 meV, A/2 = 11.18 meV, C = 0.0061 
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Figure 9. (a) Fano line shapes as a function of the Fano parameter, qp- (b) Fano 
fit to the data on (001) CeCoIns at 2.43 K. Fitting parameters used are q-p = —2.14, 
Eq = 2.23 meV, A/2 = 11.13 meV, C = 0.0061 n-\ and Go = 0.164 n-\ After [TUUj . 

Q"^, and Gq = 0.164 Q^^. In our previous paper [75], we already reported a Fano 
resonance fit to the normalized conductance data over a narrower bias region using the 
following parameter values: gp = —1-9, Eq = 0.55 meV, A/2 = 2.6 meV, C = 0.03646, 
and Go = 0.8458. From these two fits, it is clearly seen that the parameter values vary 
widely depending on the voltage range for fitting. Therefore, it is highly desirable to 
obtain and fit conductance data over a wide enough voltage range. 

One of the unresolved issues in our previous fit using a Lorentzian DOS [75] was 
why the resonance energy is located below the Fermi level. Theoretically, the Kondo 
resonance in Ce-based Kondo lattice systems is known to occur at or above the Fermi 
level, which is also observed in photoemission measurements on several heavy-fermion 
systems |134j (cf. see [129j for the Kondo resonance energy position in single impurity 
STM experiments). In the current Fano fit, it is seen that our data can be reproduced 
by the Fano formula with the resonance energy above the Fermi level, in agreement 
with theoretical predictions and other experiments. This was possible with the choice 
of a negative value for q-p, which forces the resonance energy, Eq, to be positive. In 
general, qp is a complex number since the Fano resonance arises from an interference 
effect, as discussed above. Thus, the negative qp value we obtained can be understood 
as due to the phase factor in the interference. Similar analysis using negative qp was 
also reported in some single impurity STM experiments |123j . Recently Yang [135] 
reported results from an extensive fit using a modified Fano formula to our previously 
reported data [75] over the ± 50 mV range at two temperatures (1.47 K and 20.62 K) 
and data over the ± 4 mV range at several temperatures between 2.5 K and 60 K. The 
calculated curves fit to the data well, but the implications of the extracted parameter 
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values remain to be investigated as their values depend on the range of voltage used 
in the fitting, as demonstrated above. Also, the development of a microscopic model 
beyond phenomenological arguments is needed, as was done for the single impurity STS 
study |124j . We discuss relevant issues in the following. 

One might raise questions about the validity of our analysis based on a Fano effect. 
First, how is it possible to observe a Fano resonance from point-contact (not tunneling) 
measurements and from the Kondo lattice systems? Regarding this issue, we note that 
there is an example |131j where continuous evolution of the Fano resonance across a 
tunneling-metallic junction boundary is observed from STM measurements on single 
adatoms. This implies that it may also be possible in PCS. However, one should note 
a configurational difference between the Kondo lattice and single adatoms since in the 
latter the Fano (Kondo) resonance is widely known to be sensitive to the environment 
around the impurity atom. Second, what is the origin of the interference? In the two- 
fluid picture of heavy fermions [76], there are two components in the electronic spectra 
in CeCoIns: one is the heavy electron liquid and the other is conduction electrons that 
do not participate in the hybridization. Then, one can imagine that the interference 
occurs between these two channels [751 1135] . According to the theoretical derivation 
by Madhavan et al. [124j . the Fano parameter gp is given by the ratio of two matrix 
elements: gp = ^/-B, where A is the coupling to a local orbital either direct or indirect 
via hybridization and B is the coupling to the conduction electron continuum. It is to be 
investigated further which of the two terms for A is dominant in PCS of the Kondo lattice 
system. As seen in the above fit, there exists a large background term. Go, which might 
indicate that there exists substantial contribution from the channel into the conduction 
band without causing the interference. The Fano line shape in our data persists over 
a wide temperature range including both normal and superconducting states, implying 
that the interference must be persistent over a wide temperature range. Third, how can 
we combine a Fano resonance with AR or tunneling in heavy-fermion superconductors? 
As seen in figure [5], the background conductance asymmetry affects the shape of the AR 
conductance. Thus, a successful microscopic theory should also explain this connection. 
In regards to this issue, we note that Flint et al. recently proposed a theory on the 
AR into a composite paired superconductor [136j. They claimed that it is possible 
for electrons to cotunnel into the Kondo lattice and this would provide an enhanced 
contribution to AR with a Fano resonant structure. It is interesting to look into the 
similarities and differences between their Fano resonant structure and the asymmetric 
conductance shape in our model. 

10. Remaining issues and future work 

Boundary conditions in the BTK problem. Although the theory by Deutscher and 
Nozieres [12| explains how AR is possible in heavy fermions, what role is played by the 
Fermi velocity mismatch is not clear. Related to this topic, Araiijo and Sacramento 
[112] reported recently that the argument by Deutscher and Nozieres can be proved 
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explicitly by solving this problem for a normal-metal/heavy-fermion-superconductor 
interface. Quantitative measure of the Fermi velocity mismatch effect in PCS might 
be achieved by measuring AR conductance from a ballistic junction involving a pair 
with well-matched Fermi velocities, e.g., Au/Nb, but without any oxide layer formed at 
the interface. 

Tunneling cone effect. In order to detect ABS unambiguously, making 
measurements in the tunneling limit will provide more flexible diagnostics. Thin film 
growth of CeCoIns has been attempted but the quality is not yet good enough for all- 
thin-film based tunnel junctions. As an alternative, it is worthwhile trying to deposit 
an artificial tunnel barrier onto single crystals. 

AR in heavy-fermion-metal/ conventional superconductor. Despite some success 
reported in CeCoIns/Nb point contacts, the magnitude of AR signal in such junctions 
is not well determined due to the intervening oxide layer on the Nb tip ^137j . 
Further measurements on such junctions but without oxide interface layers will provide 
important clues to many issues on AR in heavy fermions. A related experiment will 
be exploring whether proximity effect exists in such junctions |138[ 1139] . The fact that 
AR is frequently observed in heavy fermions implies that proximity effect should exist 
because it is well known that AR is at the heart of this phenomenon [140j . 

Fano effect in Kondo lattice systems. Although our proposed model has shown 
some evidence for this effect, it is not yet studied theoretically whether such effect could 
exist and, if so, is measurable by PCS. Also, it is an open question how to apply this 
picture to explain the reduced AR as well as the conductance asymmetry. Microscopic 
derivation of this effect may lead to further insightful investigations of the Kondo lattice 
physics. 

11. Conclusions 

In conclusion, PCS on CeCoIns has not only provided spectroscopic information on the 
superconducting order parameter but also clarified several experimental and theoretical 
issues. Sets of reproducible conductance spectra were obtained as a function of 
temperature, magnetic field and crystallographic orientation. A qualitative analysis of 
the data based on the extended BTK model shows that the order parameter in CeCoIus 
has a dx2-y^ symmetry. Detailed analysis of the conductance asymmetry led us to a 
phenomenological model based on a possible Fano effect in this Kondo lattice system. 
Further experimental and theoretical investigations of this model will lead us to more 
deepened understanding of the Kondo lattice physics. 
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